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ABSTRACT 
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The  Superposition  Principle  Applied  to  the  Prediction  of 
Ship  Behavior 

1.H1  — — —  ■  | 


An  outline  is  given  of  the  procedure  usually  followed  in  pre¬ 
dicting  the  response  of  a  ship  to  an  irregular  storm  sea  on  the  basis 
of  the  superposition  of  its  responses  to  all  the  components  which  make 
up  that  sea#  A  non-dimensional  form  of  presentation  is  then  developed 
which  makes  it  possible  to  visualize  and  interpret  the  results  of 
comparative  calculations  for  different  ships. 


4/0 


Acceleration  due  to  Vertical  Motions  of  a  Shi 


The  vertical  acceleration  of  a  point  in  a  ship  subject  to 
pitching  and  heaving  motions  is  considered  by  combining  vectorially  the 
pitching  and  heaving  accelerations.  A  method  of  estimating  the  expected 
extreme  values  in  irregular  seas  is  then  presented  and  compared  with 
full-scale  ship  data. 
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THE  SUPERPOSITION  PRINCIPLE  APPLIED  TO  THE 


PREDICTION  OF  SHIP  BEHAVIOR  * 

By 

Edward  V.  Lewis 

Introduction 

It  has  been  a  number  of  years  since  St.  Denis  and  Pierson 
demonstrated  the  application  of  the  superposition  principle  to  the  study 
of  ship  behavior  in  irregular  seas.  The  original  paper  (1)**  was  some¬ 
what  difficult,  but  subsequent  work  has  clarified  the  application  of  the 
technique  (2),  verified  its  validity  when  applied  to  model  and  full-scale 
data  (3)  (4),  and  developed  new  aspects  (5).  It  has  been  shown  that  the 
difficult  case  of  the  short-crested  irregular  sea  can  be  conveniently 
treated  in  a  simpler  manner  than  originally  proposed  (6). 

A  typical  plot  is  given  in  Fig.  1,  showing  how  an  oceanogra¬ 
pher's  sea  spectrum  plotted  on  the  basis  of  frequency,  CJ  ,  can  be  con¬ 
verted  to  frequency  of  encounter,  <Ue  ,  and  then  multiplied  by  a  "response 
amplitude  operator"  (for  pitch,  bending  moment,  etc.)  to  obtain  a  response 
spectrum,  assuming  a  long-crested  sea.  The  response  spectrum,  of  course, 
provides  the  basis  for  determining  various  statistical  properties  of  the 
ship  response.  The  reason  for  working  on  the  basis  of  frequency  of  en¬ 
counter  has  been  that  model  tests  in  regular  waves  show  the  peaks  of  the 
response  curves  (pitch  or  heave)  occurring  in  the  vicinity  of  the  ship's 
natural  frequency  of  oscillation. 

If  a  short-crested  sea  is  considered,  a  family  of  wave  spectra 
is  needed  to  indicate  the  energy  in  different  angular  bands  of  wave  com¬ 
ponents;  there  will  be  a  separate  response  amplitude  operator  curve  for 

*  Presented  to  the  Seakeeping  Session,  10th  International  Towing  Tank 
Conference,  London,  September  1963. 

**  See  references  listed  at  the  end  of  the  paper. 
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each  incremental  angle;  and  there  will  be  a  family  of  response  spectra. 
The  latter  can  be  integrated  over  direction  to  obtain  a  single  response 
spectrum  whose  area  and  shape  again  define  the  statistical  character¬ 
istics  of  the  .response  (6). 

As  the  results  of  systematic  model  tests  have  become  available 
(7),  it  is  now  for  the  first  time  possible  to  predict  trends  of  ship 
behavior  in  realistic  sea  conditions  (8).  However,  it  has  also  become 
clear  that  the  oceanographers'  way  of  plotting  spectra  as  a  function  of 
frequency  makes  it  difficult  to  understand  and  interpret  the  results  of 
the  calculations.  In  comparing  the  behavior  of  ships  of  different  hull 
form,  proportions,  sizes,  and  speeds  in  irregular  seas,  the  present 
system  submerges  all  of  these  factors  in  such  a  way  that  their  individual 
effects  cannot  be  separated  and  significant  comparisons  cannot  be  made. 

For  example,  Fig.  1  shows  the  comparison  of  a  250  ft.  and  a 
500  ft.  ship,  geometrically  similar  and  operating  at  the  same  Froude 
Number  in  the  same  seaway.  This  form  of  plotting  involves  first  of  all 
a  difference  in  the  sea  spectra  for  each  ship  when  plotted  on  the  basis 
of  frequency  of  encounter,  *4 •  This  is  because  the  two  ship  speeds  are 
not  the  same  at  constant  Froude  Number.  (If,  on  the  other  hand,  the 
ships  were  compared  at  the  same  absolute  speed,  then  completely  unrelated 

response  operators  would  have  to  be  used) . 

The  figure  also  shows  that  the  response  operators  have  a  very 
different  appearance,  even  when  the  ships  are  geometrically  similar  and 
the  Froude  Numbers  are  the  same.  Hence,  although  the  final  comparison 
of  response  spectra  is  significant,  there  is  no  way  to  judge  the  results 
or  to  account  for  the  differences  displayed. 

Furthermore,  a  simple  change  of  ship  speed  for  one  ship 
requires  a  transformation  of  the  sea  spectrum  and  makes  it  difficult  to 
predict  the  result.  An  obvious  modification  is  to  plot  on  the  basis  of 
(9).  In  effect  this  means  the  elimination  of  the  influence  of  ship 
speed  on  frequency,  and  the  use  of  a  frequency  scale  then  loses  its 
significance.  A  completely  new  form  of  presentation  is  clearly  needed, 
oriented  to  the  needs  of  the  designer. 
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Non-dimensional  Representation 


In  order  to  offer  a  solution  to  the  problem  of  presenting  data 
in  a  form  suitable  for  the  use  of  ship  designers,  a  non-dimensional 
scheme  will  be  presented  step  by  step  and  some  interesting  conclusions 
drawn  therefrom.  Strictly  speaking  the  proposed  scheme  is  not  com¬ 
pletely  non-dimensional.  The  term  is  used  here  only  in  the  sense  that 
geometrically  similar  ships  will  have  similar  responses  to  an  ideal  sea¬ 
way  having  the  same  component  slopes  at  all  frequencies. 

First  of  allj  we  may  recognize  that  the  parameters  describing 
ship  performance  should  be  non-dimensional.  Pitch  is  a  satisfactory 
measure  of  angular  motion,  for  it  will  be  the  same  for  ships  of  differ¬ 
ent  size  in  comparable  situations.  But  heave  or  vertical  bow  motion,  s, 
should  be  non-dimensionalized.  This  can  best  be  done  by  relating  these 
quantities  to  a  ship  dimension  such  as  length,  L,  giving  a  ratio  such 
as  s/L.  This  means  that  we  consider  two  ships  to  have  equivalent  heav¬ 
ing  behavior  in  comparable  conditions  if  the  ratios  of  heave  amplitude 
to  ship  length  are  the  same.  (This  is  in  contrast  to  the  conventional 
procedure  of  comparing  heave  amplitudes  directly). 

Once  we  adopt  a  non-dimensional  response  parameter  (an  angle 
or  a  ratio)  the  most  significant  wave  characteristic  to  use  is  the  wave 
slope,  which  is  also  non-dimensional.  This  means  that  for  the  present 
purpose  a  wave  slope  spectrum  would  be  more  significant  than  an  ampli¬ 
tude  spectrum.  Such  a  slope  spectrum  would  show  the  relative  amplitude^/ 
length  ratios  of  component  waves,  and  it  will  be  shown  that  it  can 
assist  greatly  in  visualizing  and  comparing  ship  behavior  in  an 
irregular  sea.  Strictly  speaking  maximum  wave  slope  is  A  » 

where  /.  is  wave  amplitude,  -3/'.  is  wave  height,  and  A  is  wave 
length.  Here  the  term  wave  slope  is  intended  to  mean  the  ratio  /e//A . 

The  question  then  arises  as  to  what  basis  for  plotting  sea 
spectra  and  ship  response  would  be  more  suitable  than  wave  frequency  or 
frequency  of  encounter?  Consideration  of  the  behavior  of  ship  models  in 
regular  waves  leads  to  the  realization  that  the  wave  length,  A  ,  in  re¬ 
lation  to  ship  length,  L,  is  of  equal  importance  to  frequency  of  encoun¬ 
ter  in  determining  a  model's  response  to  waves.  Furthermore,  one  could 
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make  a  non-dimensional  plot  of  response  on  the  basis  of  the  ratio. 

Such  a  diagram  would  apply  to  any  size  of  ship  or  model,  provided  the 
Froude  Number  was  constant.  However',  for  predicting  ship  behavior  in 
actual  storm  seas,  which  are  absolute  in  size,  we  must  have  a  scale  which 
directly  represents  wave  frequency  or  length,  not  a  ratio. 

A  little  experimentation  shows  that  a  suitable  scale  is  loge  A. 

If  at  the  same  time  we  use  a  non-dimensional  parameter  for  ship  response 
operator,  such  as  pitch  angle/wave  slope,  we  find  that  the  response 
operator  curves  for  geometrically  similar  ships  of  different  size  (at  the 
same  Froude  Number)  will  be  identical  in  shape,  but  separated  by  an 
amount  equal  to  log,  4  A  (Fig-  3  (b)  Furthermore,  points  at  cor¬ 
responding  values  of  ^/c  will  have  the  same  ordinates.* 

The  reason  that  a  log  scale  makes  the  response  operators 
identical  for  different  sizes  of  ships  is  that  it  makes  the  increments 
cf  ( ioge. A)  the  same  in  terms  of  ^/i.  f or  different  sizes.  For 
example,  the  increment  of  log.  A  from  1.0L  to  1.1L  is  the  same  dis¬ 
tance  on  the  log  scale  for  a  250  ft.  ship  as  for  a  500  ft.  ship.  (The 
distance  is  log-l.lL  -  log.l.OL  -  log^l.l).  It  will  be  shown  that  this 
is  a  very  desirable  characteristic  in  dealing  with  wave  spectra. 

The  relationship  between  this  proposed  system  and  the  old  can 
now  be  derived.  Since  A  -  %/■  ,  where  c  -  -  202  (ft.,  sec. 

system),  we  have! 

Cog*  A  -  loge  C  -  2( loge  CO) 

6(loge  X)  -  x  </&>--  -§-^»  °<  c 

Hence,  when  we  make  6(U constant  in  the  new  system,  the  corres¬ 
ponding  increments  in  the  old  system  vary  with  frequency.. 

Typical  response  amplitude  operators  on  the  proposed  basis  are  shown  in 

*When  this  paper  was  presented  at  the  10th  International  Towing  Tank 
Conference,  London,  September  1963,  it  was  pointed  out  by  Dr.  Y.  Yamanouchi 
of  the  Ship  Research  Institute,  Tokyo,  that  the  plot  could  also  be  based 
on  log.  instead  of  loge  A.  .  This  has  a  precedent  in  other  analogous 

work  in  the  field  of  communication  theory  and  might  be  preferred  by  some. 
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Fig.  3  (b)  for  pitching  motion.  Note  that,  although  the  curves  for 
different  sizes  of  ships  are  geometrically  similar,  they  both  fair  into 
a  value  of  1.0  and  extend  to  infinity.  This  means  simply  that  in  very 
long  waves  maximum  pitch  angle  is  equal  to  maximum  wave  slope. 

The  next  question  is  what  form  will  the  wave  spectrum  take  when 
it  is  plotted  on  the  basis  of  logfi  X  instead  of  frequency.  The  answer 
is  that  the  new  spectrum  is  related  to  the  conventional 

energy  spectrum  (in  Pierson’s  notation)  thus: 

[ 7%  >■)]*  ■  -  £ 

See  Fig.  2  (b) .  (The  minus  sign  merely  indicates  a  reversal  of  positive 
direction^)  To  prove  the  above  relationship  it  is  simply  necessary  to 
establish  that  differential  areas  (representing  squares  of  component 
wave  amplitudes)  are  the  same  in  the  two  systems.  Writing  a  differential 
area  in  the  new  system,  we  have, 

[V/oa  y] *<!((&  4 

Making  substitutions  from  above  leads  to 

[r(lo&X)¥<f  (log e  X)  = 

-  g  (-£■  &<S)  = 

The  latter  is  a  differential  area  in  the  original  system,  and  thus  the 
requirement  is  satisfied.  (See  Figure  2).  Consequently  the  area  under 
the  two  spectra  must  be  the  same;  i.e.: 

=j[r^)]  da 

In  order  to  permit  the  direct  comparison  of  ships  of  different 
size,  this  form  of  wave  spectrum  is  still  not  satisfactory.  For  this 
purpose,  as  previously  noted,  the  important  consideration  is  the  slope, 
or  amplitude/length  ratio,  of  component  waves.  Hence,  it  is  of  interest 
to  divide  points  on  the  spectrum  by  the  appropriate  wave  length  squared, 
in  order  to  obtain  a  slope  spectrum.  See  Fig.  2  (c).  Thus  a  differential 
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represents  the  component  wave 

2 

i.e.  (amplitude/ length)  . 


slope  squared  --  except  for  a 


factor  fa*)  - 


Expressing  the  above  slope  spectrum  in  terms  of  d  we  obtain, 


In  other  words,  a  spectrum  which  is  proportional  to  /cf  when  plotted 
on  the  basis  of  eJ  ,  would  be  a  horizontal  straight  line  on  the  new  pro¬ 
posed  basis  and  would  represent  an  ideal  sea  in  which  all  wave  components 
in  terms  of  logfX  had  the  same  slope.  Or  strictly  speaking  the  mean 
square  slopes  of  all  the  components  in  any  increment  of  loge*-  would  be 
the  same.  It  is  of  interest  to  note  that  there  is  some  evidence  that  the 
wave  spectrum  should  be  proportional  to  fef  under  the  influ¬ 

ence  of  a  strong  steady  wind  —  particularly  at  the  high  frequency  end 
(10).  Actual  sea  spectra  in  strong  winds  are ‘found  generally  to  be  less 
steep  than  this  (11),  as  well  as  to  fall  off  rapidly  at  low  frequencies 
(long  wave  components).  (See  the  example  plotted  in  Fig.  2).  But  the 
spectrum  proportional  to  /u)S  is  particular  interest  as  an  upper 
limit  of  sea  spectra,  in'  which  we  can  determine  how  geometrically  similar 
ships  behave.  Note  that  the  sharp  peak  in  Fig. 2a  disappears  in  2c  (arrows). 

Finally,  we  may  multiply  the  wave  slope  spectrum  (Fig.  3a)  by 
the  pitch  response  operators  (Fig.  3b)  to  give  the  non-dimensional 
response  spectra  (Fig.  3c).  These  non-dimensional  response  spectra  are 
of  direct  quantitative  significance,  since  they  represent  (pitch  ampli¬ 
tude)2,  and  the  mean  pitch  amplitude  will  be  a  function  of  the  area  under 

the  curves. 

Comparing  the  two  responses,  we  can  see  that  in  this  case  the 
shape  of  the  sea  spectrum  has  little  to  do  with  the  result.  For  even  if 
the  sea  spectrum  were  a  horizontal  line  (constant- slope)  there  would  be 
little  difference  in  response  above  logeX«  6.5.  Hence,  the  larger 
pitching  response  of  the  small  ship  results  almost  entirely  from  its 
smaller  size  --  as  indicated  by  the  horizontal  separation  of  the  response 
curves  --  which  adds  the  shaded  area  to  the  response  spectrum. 

A  similar  treatment  of  heaving  motion  leads  to  the  rather  start¬ 
ling  results  shown  in  Fig.  4.  Again  the  response  operators  are  geomet- 
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rically  similar  and  displaced  horizontally  by  a  constant  amount 
(lOgftL  /L^  ■  loge2).  But  in  this  case  they  increase  without  limit  as 
logeX  increases.  This  puts  an  even  more  severe  penalty  on  the  smaller 
ship,  no  matter  what  the  shape  of  the  sea  spectrum. 

However,  it  must  be  recognized  that  heave  amplitude  is  in 
itself  of  no  great  significance  in  ship  behavior  and  that  the  associated 
vertical  accelerations  may  be  of  much  greater  importance.  Accordingly, 
a  plot  has  been  prepared  of  non-dimensional  heave  acceleration,  a/g* 

(Fig.  5).  In  this  case  the  heave  acceleration  operators, 

mr 

must  fair  into  curves  representing  amplitudes  of  vertical  wave  accelera¬ 
tion  at  very  long  wave  lengths.  (Value  at  infinity  is  "ilP  ). 

The  resulting  response  curves  (Fig.  5)  are  surprisingly  similar 
to  those  obtained  for  pitch  (Fig.  3).  And  the  same  conclusion  can  be 
drawn:  the  larger  response  of  the  small  ship  results  almost  entirely 
from  its  smaller  size.  And,  again,  a  smaller  ship  will  always  be  worse 
off  than  a  similar  larger  ship  in  most  typical  seas. 

Usually  the  vertical  acceleration  due  to  the  combined  effect 
of  heave  and  pitch  at  a  particular  location  in  the  ship  is  of  more 
importance  than  heave  acceleration.  But  this  can  be  plotted  in  the 
same  way  and  can  be  expected  to  give  a  similar  picture  to  Fig.  5. 

Another  aspect  of  ship  motion  that  is  usually  of  even  greater 
significance  than  pitch  angle  is  the  motion  of  the  bow  relative  to  the 
water  surface  --  or  simply  "relative  bow  motion."  (See  Fig.  6).  This 
can  be  handled  in  the  same  manner  as  other  motions,  but  there  is  a  signi¬ 
ficant  difference.  As  indicated  in  Fig.  7  the  response  amplitude  operator 
for  relative  bow  motion  tends  to  zero  at  moderately  large  values  of 
log^X.  Hence,  it  can  be  seen  that  the  comparative  performance  of  ships 
of  different  size  depends  mainly  on  the  shape  of  the  sea  spectrum  and 
not  on  ship  size,  per  se. 

If  one  imagines  the  ideal  limiting  case  of  a  horizontal  slope 
spectrum,  the  behavior  of  the  two  ships  would  be  identical.  This  is  a 
case  of  truly  non-dimensional  behavior  in  irregular  seas  and  contrasts 
sharply  with  the  situations  previously  discussed  of  pitch  angle  and 
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heave  acceleration.  If  the  sea  spectrum  corresponds  to  a  decreasing  slope, 
the  large  ship  will  have  a  definite  advantage  over  the  small  ship.  In 
the  general  case,  the  comparative  performance  of  two  similar  ships  will 
depend  mainly  on  the  values  of  the  spectrum  ordinate  at  the  peaks  of  the 
two  response  operators  (approximately  at  logfcA.  ■  loggl.2  L). 

Finally,  consideration  should  be  given  to  the  response  of  a 
ship  in  terms  of  wave  induced  bending  moment.  It  is  convenient  to  use 


a  non-dimensional  coefficient, 

m  Bending  Moment 

cpg  L  3  3 

where  h  is  "effective"  wave  height,  i.e.  the  height  of  a  wave  of  ship 
length  which  would  give  the  desired  bending  moment  by  conventional  static 
calculation.  B  is  the  beam,  is  the  waterplane  coefficient,  f  is 

density,  g  is  the  acceleration  of  gravity,  and  the  coefficient  c  is  the 
static  bending  moment  coefficient  in  the  relation, 

Static  bending  mt.  ■  cfgjL  3hC„ 

Then  the  bending  moment  response  operator  may  be  expressed  in  non- 

2 

dimensional  form  as  / hg/L  \ 


From  the  sample  plot  in  Fig.  8  it  may  .be  seen  that  the  wave  spectrum  shape 
has  a  large  effect  on  the  response  spectrum  as  in  the  case  of 


relative  bow  motion. 

The  same  type  of  plotting  can  also  be  applied  for  the  case  of 
directional  sea  spectra.  Instead  of  a  single  wave  spectrum  for  each  case 
there  will  be  a  family  of  them  representing  different  component  wave 
directions.  Similarly  there  will  be  a  family  of  response  operator  curves 
and  a  family  of  response  spectra.  The  mean  square  amplitude  of  response 
will  be  a  function  of  the  double  integration  of  this  family  of  curves. 


Conclusions 

Reviewing  the  proposed  system  of  representation  for  sea  spectra 
and  the  response  o£  different  ships  thereto,  the  following  features  may 
be  noted: 
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1.  The  sea  spectrum  is  presented  in  a  manner  which 
shows  clearly  the  relative  importance  of  different 
components  in  relation  to  ship  behavior.  It  shows 
directly  the  relative  slopes  (squared)  of  the  com¬ 
ponent  waves,  which  are  more  significant  for  ship 
motion  than  wave  energy  or  amplitudes,  per  se. 

2.  The  logeA-  scale  provides  a  system  in  which  the 
response  amplitude  operators  for  geometrically  similar 
ships  are  identical  in  shape  but  shifted  along  the 
axis  in  accordance  with  absolute  size,  whereas  a 
change  in  hull  form  or  proportions  will  show  as  a 
change  in  shape  of  the  response  operator.  A  direct 
comparison  of  response  operators  for  different  ships 
plotted  in  this  manner  will  permit  a  judgment  to  be 
made  of  the  relative  merits  of  the  ships. 

3.  The  plots  of  the  final  response  spectra  give  a 
direct  visual  indication  of  the  performance  of  differ¬ 
ent  ships  or  of  the  same  ship  at  different  speeds. 

Furthermore,  one  can  see  by  a  glance  at  the  other 
curves  whether  the  superior  performance  of  one  ship 

is  the  result  of  better  motion  characteristics 
(response  operator  curves)  or  of  its  larger  size. 

It  has  been  shown  that  for  relative  bow  motions  the  ideal 
limiting  case  of  a  constant  slope  spectrum  results  in  identical  responses 
for  geometrically  similar  ships  of  any  size,  i.e.  a  truly  non-dimensional 
irregular  wave  situation.  Actual  sea  spectra  do  not  have  such  constant 
slopes,  except  perhaps  at  the  short  wave  length  (high  frequency)  end. 
Nevertheless,  plotting  of  actual  sea  spectra  on  this  basis  assists  in 
showing  the  relative  effect  on  ship  behavior  to  be  expected  from  different 
wave  length  components. 

In  the  case  of  pitching  and  heaving  it  has  been  shown  that  size 
in  itself  has  a  tremendous  effect  on  the  performance  of  geometrically 
similar  ships,  regardless  of  the  shape  of  the  sea  spectrum. 


-10- 


The  proposed  method  of  presentation  does  not  change  the  cal¬ 
culated  prediction  of  ship  behavior  in  irregular  seas  in  any  way.  But 
it  is  hoped  that  it  will  aid  in  the  understanding  and  interpretation 
of  the  superposition  technique  and  thereby  will  encourage  its  use  among 
Naval  Architects  and  research  people. 
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ACCELERATION  DUE  TO  VERTICAL  MOTIONS  OF  A  SHIP  * 


By 

Rutger  Bennet 

The  vertical  acceleration  at  a  certain  point  in  a  ship  is 
composed  of  components  from  three  modes  of  motion!  heaving,  pitching 
and  rolling.  Rolling  in  general  gives  the  least  contribution,  and  the 
greatest  part  is  caused  by  heaving  if  the  point  is  near  amidships,  and 
by  pitching  if  it  is  near  the  ends.  The  acceleration  amplitude  at  any 
point  is  very  much  dependent  on  the  phase  angle  between  heaving  and 
pitching. 

"Pitching  Axis” 

Considering  only  heave  and  pitch,  the  motion  of  a  ship  is  con¬ 
sidered  to  be  made  up  of  the  vertical  translation  (heave)  of  the  center 
of  gravity,  Zq,  and  a  rotation  (pitch, 9).  In  general  there  is  no  point 
which  during  any  finite  length  of  time  has  zero  velocity  or  acceleration; 
hence  there  is  usually  no  real  pitching  axis. 

The  motion  may  be  represented  by  a  vector  diagram  as  follows: 


\ 

o) 

/ 


The  vector  x  •  e  is  the  vertical  displacement  at  distance  x 
from  amidships  at  an  instant  when  the  angle  is  0  .  The  various  amplitudes 


*  Based  on  Report  No.  22  of  the  Swedish  Shipbuilding  Research  Foundation, 
"Accelerationer  pa  Grund  av  Fartygs  Rorelser  i  Sjo,"  by  Rutger  Bennet, 
1961  (in  Swedish). 
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are  related  by:  ** 

Z*  ■  Z*  *  * g,  X.  e  cos  e,  (D 

where  ■  amplitude  of  resulting  motion  at  distance  x  from  C.G. 
Ze  -  heave  amplitude  at  C.G.  (  X.  ■  O  ) 

Q  -  amplitude  of'.angular  motion 
£0  ■  phase  lag,  heave  after  pitch  Q 
The  minimum  amplitude  of  vertical  motion, 


A* 


Z  S/fJ  €a 
0  0 

Q 


Szebehely  (Ref.  1)  calls  this  point  of  least  vertical  motion  the  "apparent 
pitching  axis."  For  -  <  6„  <  %  ,  X  is  ne8ative> 

meaning  that  the  apparent  pitching  axis  is  aft  of  CG,  which  is  usually 
the  case. 

If  both  motions  are  harmonic,  the  variables  in  equations  (1) 
and  (2)  may  be  directly  replaced  by  their  first  or  second  derivatives. 
Therefore,  the  point  with  the  least  vertical  motion  will  also  have 
the  least  vertical  velocity  and  the  least  vertical  acceleration. 

Irregular  Motions 

In  irregular  waves,  when  motions  are  no  longer  harmonic,  con¬ 
ditions  become  more  complicated.  The  point  with  least  acceleration  will 
not  in  general  have  the  least  displacement,  and  none  of  the  minimum 
points  will  be  fixed  but  will  move  between  plus  and  minus  infinity.  The 
distance  will  vary  according  to  a  statistical  distribution  with 

mean  value  according  to  equation  (2)  (see  Ref.  2).  It  is  now  generally 


**  See  summary  of  Notation  at  the  end  of  this  report. 
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accepted  that  all  ship  motions  may  be  represented  by  stochastic  processes 
which  are  Gaussian  or  normal.  This  means  that  the  distribution  of  equally 
spaced  values  of  Z  and  9  should  be  normal;  thus: 


(3) 


with  a  similar  expression  for  yO(&). 

The  amplitudes  should  be  Rayleigh  distributed: 


a. 


6) 


where  a  is  the  amplitude  of  any  motion  and  R  is  the  mean  square  amplitude. 

At  any  point  in  the  length  of  the  ship,  the  sum  of  accelerations 
due  to  pitch  and  heave  will  also  be  normal,  with  Rayleigh  distributed 
total  amplitude.  The  parameters  of  the  distributions  will  be  functions 
of  the  correlation  between  pitch  and  heave  (Ref.  2). 


Energy  Spectrum 

Any  irregular  motion,  e.g.  pitching,  may  be  written  as  a  sum 
of  a  large  number  of  regular  components: 


Q(t)  -  £  **  co3(cJ„t  +  en) 


(S) 


The  energy  spectrum  is  defined  as  the  spectral  density  in  a 
frequency  interval  cT*sJ  ' 


ekj 


4  ' 


£4 

&C<) 


(6) 


For  each  component  the  acceleration  is  the  second  derivative: 

4  (t)  =  -  ® 

The  acceleration  spectrum  is 


(s) 
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The  variance  of  the  statistical  distribution  is  the  area  of 
the  spectrum;  and  this  is  also  equal  to  half  the  parameter  R  in  eq.  (4): 


JS*  =  (f 
z  « 


</eJ 


4)  4  dcd 


If  the  combined  effect  of  pitch  and  heave  on  the  motion  of  a 
point  is  known,  then  spectra  of  motion  and  of  acceleration  of  the  point 
can  be  obtained  in  the  same  manner.  However,  in  many  cases  the  pitch 
and  heave  motions  are  only  known  separately.  It  is  of  value,  therefore, 
to  develop  approximate  methods  of  estimating  accelerations  at  a  point 
from  pitch  and  heave  motion  data,  preferably  without  actually  computing 
the  spectra  of  either. 

Angular  Acceleration  Deduced  from  the  Angle 


If  every  period  of  the  motion  in  irregular  waves  were  harmonic, 
the  acceleration  amplitude  would  be 


4  If'*' 

^9  =  T *  ^0 


If  the  period  T  could  be  considered  as  constant  also  in  an  irregular  sea, 
there  would  be  a  simple  relationship  between  the  statistical  parameter  R 
(eq.  4)  for  pitch  amplitude  and  acceleration: 


ft  •(¥)  ft 


iis  is  equivalent  to  (J  ■  const,  in  equation  (9b). 

But  as  the  period  T  is  varying,  equation  (11)  will  not  hold  in 
1  irregular  sea.  It  is  necessary  to  apply  a  correction  factor  K,  so  that 
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or,  with  t  -  V A) 


03) 


ri 


<■ 


An  investigation  of  measurements  of  pitching  angle  and  angular  accelera¬ 
tions  in  Ref.  (4)  has  given  values  of  K  between  1.3  and  1.6.  In  ref. 

(5)  K  -  1.5  is  given  as  a  recommended  value.  If  the  spectrum  has  an 
equation  of  the  same  type  as  the  Nevmann  wave  spectrum,  K  -  1.5  is 
obtained  exactly. 

In  Ref.  (3)  spectra  of  pitching  angle  for  a  cargo  ship  are 
presented,  and  for  two  of  these  the  acceleration  spectra  have  been 
derived  according  to  equation  (8).  From  the  relation  between  the  areas 
of  the  angle  and  acceleration  spectra,  K  was  found  to  be  1.2  and  1.35, 
respectively. 

It  may  be  concluded. that  the  acceleration  caused  by  an  ir¬ 
regular  pitching  motion  can  be  up  to  507.  greater  than  in  a  regular 
motion  with  the  same  amplitude. 


Vertical  Acceleration  from  Heave  and  Pitch 

At  a  distance  X  from  C.G.  the  amplitude  of  vertical  accelera¬ 
tion  according  to  eq.  (1)  is! 

a,x  -  V K  +  *  *■*'•*'*•#  co5  £o 

In  terms  of  a  known  acceleration  0.^  at  a  point  Xm  A  this 


equation  becomes 


a*  -  y  ”■  ***(*--4)  «t cos 


(!*) 


(a) 


i*'  A 

where  C.  is  jJ^ase  lag  between  heave  and  pitch  at  X  “  f • 
A 

,'The  time  phase  angle  is 

a.A  5/m  £a 

£  =  4RCTA(J  — 2 - - - 

0  Ol.  Cos 


m 


In  analyzing  measurements  it  is  an  advantage  to  use  mean  values 
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instead  of  individual  amplitudes.  The  mean  square  of  the  resulting 
amplitude  is 

(t-ATK  *-<Hcor£«  (/?> 

The  last  term  is  the  mean  of  a  three-dimensional  distribution.  If  the 
three  variables  are  independent,  the  mean  of  the  product  is  the  product 
of  the  means: 


Calculations  and  measurements  have  shown  that  the  coupling  between  pitch 
and  heave  has  only  a  small  influence  (Ref.  2);  so  eq.  (18)  may  be 
assumed  to  be  valid. 

The  Rayleigh  distribution  has  a  constant  ratio  between  mean 
value  U  and  the  mean  square: 


where  y  is  any  Rayleighian  variable. 

Eq.  (17)  can  then  be  written  in  the  more  useful  form: 

If  the  angular  acceleration  is  not  known,  but  only  the  angle,  is 
calculated  from  eq.  (12)  with  K  -  1.5. 


Expected  Extreme  Values 

The  value  that  is  expected  to  be  exceeded  with  the  probability 
Q  is  determined  from  the  properties  of  the  Rayleigh  distribution: 

a-t0?)  '  (*0 

where  is  a  function  only  of  Q.  In  ref.  (6)  it  is  shown  that 
Q 
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[/.  ■  3.5  is  a  reasonable  extreme  value  for  practical  consideration. 

<2 

Theoretically  it  is  the  value  which  with  997.  probability  will  only  be 
exceeded  once  during  2000  cycles,  corresponding  to  about  four  hours  sail¬ 
ing  time.  This  value  will  be  used  below  to  predict  expected  extreme 

| 

accelerations. 

Results  of  Measurements 

Since  actual  data  on  observed  accelerations  are  scarce,  a 
summary  may  be  given  of  available  published  information. 

DTMB  Full  Scale  Tests 

Results  of  motion  and  acceleration  measurements,  carried  out 
by  David  Taylor  Model  Basin  on  the  SILVER  MARINER,  are  published  in  Ref. 

(4).  The  observations  were  made  at  various  speeds  and  headings  in  a 
fairly  moderate  sea  with  a  significant  wave  height  of  9.4  feet  and  a 
wind  speed  of  23-28  knots.  The  highest  recorded  heave  acceleration  at 
C.G.  was  0.26  g,  and  during  the  same  cycle  the  angular  acceleration  was 
0.11  rad/sec.2  with  a  phase  angle  of  40°.  This  was  also  the  highest 
recorded  value  of  Q  .  The  maximum  root  mean  square  amplitudes  were. 

tf**  r0  *  o.  //  g 

Afjx.  Q;  ~  0.049 

Japanese  Model  Tests 

Model  tests  are  published  in  Ref.  (6)  for  three  different  ship 

types: 

Tanker,  Displacement  87,500  tons,  Cg  »  .80, 

Cargo  liner,  Displacement  17,400  tons,  Cg  ■  .646, 

"Experimental  Ship",  Displacement  3,460  tons,  Cg  -  .52. 

The  tests  were  made  in  regular  waves,  and  vertical  accelerations  were 
measured  at  two  points. 

The  superposition  method  was  used  to  calculate  the  response  to 
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a  Neumann  wave  spectrum.  Root-mean-square  acceleration  amplitude  was 
determined  from  equation  (9),  and  expected  extreme  values  from  eq.  (21). 
Fig.  1  from  Ref.  (6)  gives  values  which  with  90%  probability  will  only 
be  exceeded  once  in  5000  cycles.  This  means  ■  3.3  in  (21),  which  is 
close  to  the  value  3.5  suggested  above.  The  Japanese  investigation  was 
made  in  order  to  determine  the  best  position  in  a  ship  for  a  nuclear 
reactor. 

Swedish  Tests 

The  Swedish  Shipbuilding  Research  Foundation  has  measured 
stresses  and  motions  in  a  cargo  liner  and  an  ore  carrier,  Ref.  (3). 

M/S  CANADA  is  a  465 * ,  19.5  knot  cargo  liner.  Pitch  angle  and 
vertical  acceleration  at  a  point  42  feet  aft  of  C.G.  were  measured. 

The  highest  recorded  acceleration  was  0.25g,  and  during  that  cycle  the 
pitch  angle  was  4°  with  90°  phase  lag.  The  highest  recorded  pitch  angle 
was  7.5°,  with  acceleration  0.21  g,  also  with  90°  phase  difference.  The 

period  T  was  6.5  sec.  assuming  K  -  1.5,  this  means  an  angular  accelera- 
2 

tion  of  0.184  rad/sec. 

From  eq.  (15)  and  (16)  is  found: 

Heave  acceleration  at  C.G.  CLa  -  0.32  g 

Phase  angle  at  C.G.  £.  -  41° 

2 

Angular  acceleration  <2.v“  0.184  rad/sec. 

The  distribution  of  this  acceleration  along  the  length  of  the 
ship  is  shown  in  Fig.  2,  labeled  "Max.  recorded." 

The  same  record  gave  the  following  root  mean  square  values: 

tg  =0.073  r"%ec* 

%  ~0.//4g 

At  C.G.  re  -  0.14  g) 

At  F  P  -  0.58g  )  Total  vertical  acceleration 

At  A  P  rx  -  0. 37g  ) 


The  number  of  cycles  in  the  record  was  300,  and  the  most 
probable  highest  value  is  then  Z.4  K  (  I/q  *  Z.4-J  .  Table  I 
shows  good  agreement  between  theoretical  and  measured  maxima.  The 
table  also  shows  tht  highest  value  expected  in  2000  cycles  according 
to  eq.  (21)  (V  -  3.5).  The  longitudinal  distribution  is  shown  in 
Fig.  2. 

A  limited  number  of  measurements  was  also  made  on  a  560*  ore 
carrier,  M/S  VASSIJAURE  (Ref.  3).  The  results  are  given  in  Table  II 
and  Fig.  3. 

These  examples  show  how  limited  information  on  pitching  and 
heaving  motions  will  permit  good  approximate  estimates  to  be  made  of 
maximum  expected  accelerations  at  a  particular  point. 
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TABLE  I 


M/S  "CANADA 

L  -  465',  B  -  64',  d  -  26.7', 
C_  -  0.65,  Displ.  -  14,500  ts. 

D 


r 


Theor.  max  Maximum 

ampl.  during  measured 

measurement  amplitude 


Max. expected 
amplitude  in 
2000  cycles 


Heave  acceleration  at 

C.G.,  g's 

0.14 

0.33 

0.32 

0.5 

Pitch  angular  acceler¬ 
ation,  rad/sec 

0.07 

0.175 

0.184 

0.25 

Vertical  acc.  42'  abaft 
C.G.,  g's 

0.11 

0.27 

0.25 

Vertical  acc.  at  FP,  g's 

0.58 

1.40 

1.60 

2.0 

Vertical  acc.  at  AP,  g's 

0.37 

0.90 

1.07 

1.3 

TABLE  II 

M/S 

"VASSIJAURE" 

L  -  560 

',  B  -  74.5',  d  -  31', 

r 

B  -  0.' 

78,  Displ.  -  29,000  ts. 

Maximum 

measured 

r  amplitude 

Theoretical  max. 
amplitude  during 
measurement 

Heave  acceleration  at  C.G.,  g's 

0.035  0.12 

0.07 

Pitch  angular 
rad/sec^ 

acceleration, 

0.046  0.16 

0.13 

Vertical  acc. 

at  FP,  g's 

0.42  1.45 

1.15 

Vertical  acc. 

at  AP,  g's 

0.38  1.35 

1.09 
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NOTATION 


Q  Pitching  angle 
Za  Heave  at  C.G. 

Z  -  Z„  +  X  Q  -  Vertical  displacement  at  a  distance 
x  from  C.G. 

^  x-coordinate  of  accelerometer 
Amplitude  of  pitch  angle 
Zi0  Amplitude  of  heave  at  C.G. 

Amplitude  of  heave  at  x  from  C.G. 
t L  Amplitude  of  motion  or  acceleration 
CL3  Amplitude  of  angular  acceleration 
6L6  Amplitude  of  heave  acceleration  at  C.G. 
dLy  Amplitude  of  heave  acceleration  at  x  from  C.G. 
od  Angular  velocity 

Phase  lag,  heave  after  pitch  at  C.G. 

Cp  Phase  lag,  at  x  from  C.G. 
o'  Variance  of  a  statistical  variable 
^  Mean  square  amplitude 

r  =^f7 

X.  Coordinate  in  longitudinal  direction 
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